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Using an exact mapping to disordered Coulomb gases, we introduce a novel method to study 
two dimensional Dirac fermions with quenched disorder in two dimensions which allows to treat non 
perturbative freezing phenomena. For purely random gauge disorder it is known that the exact zero 
energy eigenstate exhibits a freezing-like transition at a threshold value of disorder a = a t h — 2. 
Here we compute the dynamical exponent z which characterizes the critical behaviour of the density 
of states around zero energy, and find that it also exhibits a phase transition. Specifically, we find 
that p(E = + ie) ~ e 2/z ~ x (and p(E) ~ E 2/z ^) with z = 1 + a for a < 2 and z = V&F - f 
for a > 2. For a finite system size L < e -1 / z we find large sample to sample fluctuations with a 
typical jO e (0) ~ L z ~ 2 . Adding a scalar random potential of small variance 6, as in the corresponding 
quantum Hall system, yields a finite noncritical p(0) ~ S°" whose scaling exponent a exhibits two 
transitions, one at crth/4 and the other at o t h- These transitions are shown to be related to the one 
of a directed polymer on a Cayley tree with random signs (or complex) Boltzmann weights. Some 
observations are made for the strong disorder regime relevant to describe transport in the quantum 
Hall system. 



I. INTRODUCTION 



The critical behaviour of the plateau transitions in the 
integer quantum Hall (QH) effect remains an appealing 
theoretical challenge. Despite numerous attempts, a cal- 
culable theory remains elusive. An equivalent version 
of the quantum Hall system which is believed to cap- 
ture the relevant physics corresponds to two dimensional 
Dirac fermions in presence of both a random vector and 
a random scalar potential . Conventional perturbative 
methods have failed and it is believed that the problem 
is described by some non perturbative strong coupling 
regime QJ^]. Recent works using conformal field theory 
[^-Q or non linear sigma models aim at reaching this 
regime M. 

One possible route of attack is to use the boson repre- 
sentation JlJIJnl based on the network model In- 
deed, the model can be mapped exactly, via bosoniza- 
tion, onto a random sine Gordon model or equivalently a 
Coulomb gas (CG) with a specific type of disorder. Al- 
though the calculation of the density of states, via the 
retarded Green's function, corresponds to considering a 
single CG layer, the full treatment of the quantum Hall 
transition (both advanced and retarded Green's function) 
requires to study two coupled Coulomb gas layers and 
remains highly non trivial in these new variables. On 
the other hand, there has been recent progress in under- 
standing disordered CG, mainly in the context of ran- 
dom gauge XY models [|l2] [l7 19 and in particular the 
freezing transitions which occur in these systems. New 
methods, such as fugacity distribution rcnormalization 
group (RG) Jn|-(l8| as well as variational methods p9| , 
have been developed which seem to be able to capture 
some of the non perturbative features of the strong dis- 
order regimes. It is thus of interest to search what can 
be learned from these methods and to understand, irre- 



spective of formal technicalities, whether the (glass tran- 
sition) physics that they describe will be part of the QH 
strong disorder physics. 

In this paper we mainly focus on the detailed under- 
standing of the single layer problem in the Coulomb gas 
formulation with the practical aim of computing the den- 
sity of states. We also extend our method to the full QH 
problem, proposing a new approach on this venerable 
problem. 

We start by further restricting to the purely random 
vector potential disorder model, the scalar random po- 
tential will be added later on. This simpler model has 
been intensively studied Jl|,0,^ 20 2l]] and is believed to 
be critical, with a line a fixed points, and a continu- 
ously varying dynamical exponent z{a) as a function of 
random vector potential disorder strength a. Some pre- 
cise results exist for an exactly known zero energy eigen- 
state which has the form ip(x) = e u ^ x ^ 2 ijjo where U{x)/2 
is the primitive of the vector potential. It was found 
p2[ that averaged moments scale with system size L as 
^2 X \ip(x)\ 2q ~ L~ T ( q \ such that above a threshold value 
o = <?th of disorder r(q) = for sufficiently large q indi- 
cating some kind of localized behaviour. Further studies 
|l7f confirmed the existence of a transition at a = <r t h 
in the (Gibbs like) probability measure |^(x)| 2 = e 11 ^ 
(equivalently a freezing, i.e. a glass transition), as well as 
its relations, via RG, to the directed polymer on a Cay- 
ley tree |^3|,|l4],^0| and found a non trivial structure of the 
strong disorder phase with "quasi localized" behaviour. 
Interesting relations to the Liouville theory, conjectured 
in |2l| were reexamined and it was found that the freezing 
transition can be directly demonstrated from renormal- 
ization in the Liouville model | p7[ . 

The known results about the exact E — eigenstate 
p2[ do not however tell anything directly about the den- 
sity of states. In particular the dynamical exponent has 
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not yet been calculated in the strong disorder regime, 
and one would guess that it should exhibit some kind 
of change at the transition a = a t h- A freezing in the 
dynamical exponent was indeed demonstrated recently 
p4j in a closely related model, i.e. the classical Arrhe- 
nius diffusion in the potential U(x), in both one and two 
dimension at the same value a = ath than the E = 
eigenstate transition. In one dimension the square of the 
Dirac Hamiltonian is well known to be identical to the 
Fokker Planck operator and there the two problems are 



thus equivalent |25|. Thus in one dimension, if one con- 
siders a log correlated U(x), both problems have identi- 
cal dynamical exponents and freezing transitions given in 
Ipij . In two dimensions, as discussed below, the two dif- 
fer by a an additional imaginary random drift term, but 
still they both have a line of fixed points and it is reason- 
able that they would both undergo freezing transitions, 
as we find here. 

In this paper we start by defining the models (Section 
II) and by showing that the density of states (DOS) of 
the Dirac Hamiltonian can be expressed as an observable 
in a boson formulation. For convenience we study the 
DOS p e {E) at energy E = adding a small but finite 
imaginary ie term for the retarded propagator, thus in 
effect computing a smoothed DOS, and carefully study 
the limit e 0+ (Section III). At E = the model 
becomes very similar, in the boson formulation, to the 
random gauge XY model in the phase where vortices are 
relevant. The parameter e plays the role of a bare vortex 
fugacity and the local DOS p £ (0,r) corresponds to the 
renormalized vortex fugacity z± (r) (or the local density) 
which becomes broadly distributed when e — ► 0. We show 
that the order of limits e — > and system size L — > oo 
is significant. For L — > oo, such that the typical level 
spacing AE < e, we use a variational scheme and show 
that p e (0) ~ e 2 / z_1 with z exhibiting a transition at a 
critical value of disorder. This is equivalent to a phase 
transition in p(E) ~ E 2 / Z ~ l . For AE > e we find that 
p e (0) becomes analogous to to the partition function of 
a directed polymers on a Cayley tree, and also exhibits 
the freezing transition. It is however a strongly fluctuat- 
ing quantity in that limit and is interpreted as a typical 
value, rather than a disorder average. Further analogies 
with freezing of dynamical exponents in Arrhenius dy- 
namics is presented. Finally, in Section IV we include a 
scalar random potential as in the full Quantum Hall sys- 
tem. We find that the DOS is non critical in E, however 
its a dependence is critical. We also develop a varia- 
tional scheme for studying the transport and localization 
exponents. 



II. SINGLE LAYER MODEL, DEFINITIONS, 
AND EXACT MAPPINGS 



Our aim is first to study the density of states of the 
random Dirac Hamiltonian in two space dimensions: 



H D = hvper ■ [-iV - A(r)] + W{r) 



(1) 



where r = (x, y) is the 2D space, a = (a x ,a y ) are Pauli 
matrices, W(r) is a random scalar potential and A(r) is 
a random vector potential (in units of e/h), both gaus- 
sian with short range correlations (in the following we set 
hvp = 1). A can be chosen purely transverse A x = d y V, 
A y = —d x V and its potential has logarithmic correla- 
tions (V(r) — V^r')) 2 ~ a In |r — r'|, which defines a. 
Two exact zero energy (unnormalized) eigenstates are 
then ip = (e v , 0) and ij) = (0, e~ v ). 

Since we are also interested in the local density of states 
in a given sample and that this is a fluctuating quantity 
it is convenient to define the smoothed local density of 
states p t (E, r) as: 



p e (E,r) = -3<r|— — 

7r E — Hd — ie 



r > 



(E - E n y + e 



#n(r)| 2 



(2) 



where E n are energy eigenvalues of Hp and i\) n the eigen- 
states. For e — > + one recovers the standard local DOS, 
and for finite e each level is broaden by a Lorentzian. The 
standard DOS is then defined as the spatial average, for 
a system of linear size L: 



Pe{E) = ^ J d 2 r Pe (E,r) 



(3) 



Although this usually becomes a smooth function of E 
for L — * oo in any given sample, at finite size and for 
small e it is a series of peaks whose locations usually 
fluctuate strongly from sample to sample. Clearly these 
fluctuations are smoothed when e becomes of the order 
or larger than the typical level spacing AE. Naively, if 
the lowest energy states scale as L~ z then dimensional 
argument gives for the space averaged DOS p £ (0) ~ L z ~ 2 
for e < L~ z or p e (0) — e 2 / z_1 for e > L~ z . 

The local DOS can be expressed from the free fermion 
action with spinors ip(v), iP(t), projected into a subspace 
of energy E which defines the Dirac problem in 2 + I 
dimensions. 



f 



p e (E,r) = -9 < ^(r)ip(r) > Sd 



d 2 riP(r){(T ■ [-iV-A(r) 



+W(r) - E + ie}ip(r) 



(4) 



(5) 



An additional Dirac mass term Ad<J z in Eq. (g) controls 
the distance from criticality and is set here to zero. 
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The problem can be mapped onto a sine Gordon model. 
Considering y as an imaginary time variable this action 
can be written as a 1+1 dimensional fermion problem. 
Further bosonization jj| yields the action 

d 2 r{ ' 



S 



na 



(r)] 2 

-(W(r)-E + ie) cos6>(r)} 



+ —[A v {v)d x -A x {v)d y ]e{v) 



(6) 



where a (which denotes Hvfck) is the momentum cutoff, 
and K = 1; K ^ 1 may be generated by RG or corre- 
spond to ID type interactions. Allowing for K ^ 1 is 
mainly instructive as it allows to interpolate towards the 
random gauge XY model, and we call that situation the 
generalized Dirac model. Eq. (^J) can also be derived 
from the network model of the QHE P-pl| where x is 
discretized, W(x,t) is then the long range component of 
the random potential while A y (x, r) = ( — ) x W(x, r) is 
the short wavelength component; both terms couple to 
slowly varying fields, hence can be considered as inde- 
pendent random variables, though with equal averages. 

In the sine Gordon formulation, the local (smoothed) 
DOS is given exactly by an average of the operator cos(6*) 
as follows: 

Pe(E,r) = -^-£--mnZ=-L?(l{cos9(r)) (7) 
irL^ dW(r) ■n i a 

where Z = J Dip Difie~ Sd and (cos9) is an average over 
8 with the action Eq. (|6|) . 

We note that all the above mappings are exact. They 
are even exact for a finite size sample with some specified 
boundary conditions for the path integrals. We will not 
need here to detail the correspondences in boundary con- 
ditions, but it can be done in principle. Note that since 
the action is complex, the (cos9(r)) can be arbitrarily 
large (when the denominator vanishes) which is the case 
for e — > + as E crosses an eigenvalue E n . 



III. RANDOM VECTOR POTENTIAL MODEL 

We now set W(r) — and study the model with only a 
random vector potential. To determine the dynamical ex- 
ponent z we will study the smoothed DOS at zero energy 
E = 0. Below we will distinguish two limits and study 
them separately. First in the large size limit (e > L~ z ), 
if we assume, as is customary, that there is a well defined 
density of states p(E) ~ _E 2 / 2-1 , when L — > +oo one has: 



dEp(E)- 



2/z-l 



E 2 



(8) 



for fixed small e and z > 1. Thus we can obtain z un- 
ambiguously from p e (0). This observable should be self 
averaging in that limit since the DOS at zero energy re- 
ceives contributions from many energy levels in a window 
of size e around E — 0, and this is what we find below. 



There is another interesting limit, also studied below, 
when e < is small (respectively finite size). Then 
there are fewer energy levels and p e (0) becomes a strongly 
fluctuating quantity, as discussed below, which gives in- 
formation about the statistics of the lowest energy levels 
Eq near E = 0, and thus also about the typical energy 
level spacing, found to scale also as Eq ~ L~ z with the 
same exponent z. 



A. Large size limit 

A look at Eq. (|) setting E = (with W(r) = 0) shows 
that the model is identical to the sine Gordon (or equiv- 
alently Coulomb gas) formulation |^6|,|9| of the random 
gauge XY model. Clearly e plays the role of the bare 
vortex fugacity, but as shown in Jl6| the random vector 
potential generates, upon coarse graining, an additional 
local random potential resulting in a random fugacity 
z±(r) = ee ±c/ ( r ) for ±1 charges. Some of the physics 
of the random gauge XY model will thus be relevant 
here. In particular the local DOS p e (0,r) is analougous 
to the coarse grained z± (r) and thus become broadly dis- 
tributed as e — > 0, as discussed below. 

It is convenient to perform a replica average on A in 
Eq. (^|). This yields a Hamiltonian for the replicated 
field a (r) with replica indices a,b = 1, m 

H = J d 2 r{{l/^)Y^{K- x 8 ab + (j)Ve a -V6 b (9) 

ab 

-^Y>]exp(m-0)}. (10) 

n 

where (A 2 , (q)) = (A 2 (q)) — it<t, n is a vector of length m 

with entries 0, ±1 and Y[n] ~ J\ a e™». The term ecos(9) 
in Eq. (Jfy corresponds to Y[n] with ^ Q r? a = 1 while all 
other n are generated by renormalization group (RG). 
The inclusion of all these terms is essential for treating 
properly the strong disorder situation |T^,[l6|,[l9| , and ob- 
taining the correct scaling dimension of the e cos(#) op- 
erator, which is what we need here. Since e is finite, and 
we are mostly interested in the region K > 1/2 where 
the vortices are relevant, there will exist in finite density 
(separated by a scale e~ z ). Since we are interested in the 
end in the behaviour as e ^ (dilute limit) we can use 
the RG method developped in |16|] and follow the full dis- 
tribution of fugacities or equivalently all the Y[n], up to 
the length scale at which the vortices separation becomes 

— 1/2 

of order one, corresponding to L ~ a c below (see Ref 
p8[ for a similar RG study). 

For simplicity we use instead a variational method, 
shown in our previous studies fTsj ] to be good enough 
to describe the dilute vortex system. In the limits of 
interest in Sections III and IV A (small E, e) this varia- 
tional method is easily seen (by comparison to the above 
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mentioned RG) to give the exact result for scaling dimen- 
sions. Very much as in Jl6| we expect the more precise 
RG treatment to correct only amplitudes at weak disor- 
der, powers of logarithmic prefactors at strong disorder 
and be necessary mostly for detailed descriptions very 
near the transitions, which we leave for future publica- 
tion. The results given below for the exponents z and z' 
should thus be considered as exact. 
The variational Ho has the form 

Ho = J dM(l/8vr) ^[OFT- 1 ^ + a)V9 a ■ V6 b (11) 



ab 



(12) 



with <7 C , Co variational mass parameters. The propaga- 
tors of the 9 field are used to define 

$Z<0a(q)0b(-q)>O = - 2u5 ab - A 

q 

u = -(K/2) ln(A c / A-kK<t c ) 

A = aK 2 m(A c /47riTo- c ) + Ka /a c - aK 2 . (13) 
The interaction term is (e/na — > e here) 

(^y[n]exp(zn.0)) (14) 

n 

= (5^eaqp[(« + Inc)5^n2 + w5^n B ]) w (15) 

ii a a 

where the w average reproduces the required form with 



(16) 



. exp[-uj 2 /2A}duj/(2Trv / A) 



and has the physical interpretation of an average over 
random local fugacities [jl9|,[l6| • The sum in Eq. ([H]) can 
be written as (Hi n t)o = (Z m )u with 



Z = 1 + ee u+w + ee u ~" 



(17) 



The variational free energy is then minimized, F, 
Fo + (H — Ho)o, where F is the free energy of Eq. (|l 
and ()o is an average with weights exp(— Ho). This pro- 
cedure yields JTsf] an equation for o~ c 



= dui 



ce 



U+U) I „U — UJ 



2Jlu 



(1 + ee u+UJ + ee"-^) 5 



>,2A_fa (18) 

V2^4 



where 



u = {K/2) ln(47rifcr c /A c ) 

A = -aK 2 ln(4nKa c /A c ) - K 2 a + Ka /a c 



(19) 



and A c 3> <J C is an integration cutoff. Eq. (jlq) can be 
solved by steepest descent when the logarithms are large. 
A similar equation for ctq yields that ero/er c in Eq. (J19|) 



is at most finite and therefore can be neglected to deter- 
mine exponents. The result is a phase diagram shown in 
Fig. 1 with a massive phase <j c ^ bounded by the lines 
2 — K + aK 2 = and a = 1/8. Furthermore, the line 
a = 2/K 2 manifests a phase transition corresponding to 
a change in the relation a c ~ t 2 ^ z where 



z = 2 - K + a K 
z = K(V8a - 1) 



a < 2/K 2 
a > 2/K 2 . 



(20) 



This transition occurs as rare regions of the sample rather 
than typical ones start dominating the behavior p6]| , as 
can be also seen from ©. For a < 2/ K 2 one can discard 
denominators (as well as the e 2 term) which immediately 
yields (|20|). For a > 2/K 2 the average over the random 
fugacity u> is dominated by the tail of the distribution 
and the right hand side of (|l8|) can be approximated by 
& 



Proba(w + u + In e > 0) w exp( 



8aK 2 



lna c ) (21) 



yielding the strong disorder form of (pOj). 

The above results now allow to compute straightfor- 
wardly the disorder averaged DOS. Indeed we can identify 
the disorder average of Eq. (]|) as p e (0) — -^F var /nL 2 
with the overline denoting disorder average. This yields 
the replica average 



ld_ 

it de 



^(°) = T-(E r [ n ] ex p( m - )>- 



(22) 



Using E q. (|l8| ) for a c and a corresponding equation for ao 
we find p e (Q) - (a c + a )/e ~ e 2 /^ 1 . Note that Eq. @ 
has beyond the (cos 9) term all the higher order terms 
in e as generated by RG. The result differs from just the 
(cos0) average in the strong disorder regime, a > 2/K 2 . 

In this derivation we have used that L is large com- 
pared with the correlation length l/«/o^ ~ e -1 / 2 so that 
integration cutoffs are determined by a c . Hence our re- 
sult is that 



Pe(P) 



e > L~ 



(23) 



exhibits a phase transition at a = 2/K 2 . From (|§) this 
also implies that the DOS also exhibits a transition, with 
p(E) ~ E 2 /*- 1 . The condition e > L~ z can be inter- 
preted as the typical level spacing AE = l/p e (0)L 2 be- 
ing small, AE < e, so that levels overlap and the DOS is 
smooth at finite L. 

Our conclusion for the free random Dirac Hamiltonian 
(|lf), obtained by setting K = 1, is that there is a phase 
transition at a = 2. We find: 

p(E) ~ £(i-°-)/(i+<7) for a < 2 (24) 
p(E) ~ £j(3-V8?)/(v^-i) for a > 2 (25) 

The value obtained for the threshold coincide with the 
one for the exact zero energy eigenstate 20,22]. Thus at 
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this value of disorder a sharp change of behaviour also 
occurs in the DOS. A more detailed treatment reveals 
that in the strong disorder phase there are logarithmic 
prefactors to the DOS, as was also the case in For 
(j > 2 one has: 



p{E) ~ E 2 l z ~ 1 \\nE\ 



-4, 



(26) 



'8cr-l 



with from 



where we find for E = ie that ip — 
©7=1/2. 

As mentioned above, using the RG of [|l6| yields the 
same result (20) for the dynamical exponent z. In fine, 



this is a result about the true scaling dimension z of the 
ecos(6 l ) operator (noted A typ in the conclusion of 16 1) 
being different, in the strong disorder regime, from the 
naive one (noted A there), as occurs in the XY model 
(see discussion there). The RG treatment is expected to 
change the value of the exponent ip of the logarithmic 
corrections (i.e. 7 is expected to change to 7 = 3/2 [0). 



O 

2 



1/8 





/ 




/ 

/ 


/ 

/ 

Strong disorder 

✓ 

✓ 








✓ 

s 

s 

^1 Weak disorder 




Massless \ 





1/2 



1 



1/K 



FIG. 1. Schematic phase diagram for the generalized 
Dirac problem (e.g. single layer disordered Coulomb gas). 
a is the strength of the random gauge disorder. K is an inter- 
action parameter so that K = 1 corresponds to free fermions. 
The full line is the phase boundary above which single charges 
become relevant and below which the model is massless. The 
dashed line indicate the freezing transition between weak and 
strong disorder regimes. 



B. Finite size regime 

Let us now characterize some aspects of the fluctua- 
tions of the DOS in a finite size system. Within the vari- 
ational method described above one sees that for e > L~ z 
the system is too small to generate the mass a c hence 
/9 e (0) ~ (cos 0)o with ct c = (T = in Eq. (O), i.e. 



Pe(0) 



-K+oK A 



e < L~ 



(27) 



Since the system is effectively massless we expect signifi- 
cant fluctuations. In the following we consider a different 



approach for the e < L~ z case which will clarify the na- 
ture of disorder average. 

We proceed to evaluate the DOS by a direct expansion 
in e. At e — a direct evaluation of the Gaussian aver- 
age over in a given sample (assuming periodic boundary 
conditions for the resulting potential V(r)) yields [p2[: 



(cos8(r)) 



e -KlnL le -U(r) 



,U{r) 



]/2 



(28) 



where, in Fourier space, U{q) = 2KV(q) = 
(2K/q 2 )(iq x A y — iq y A x ) with correlation of the form 



(C/(r) - C/(0)) 2 = AaK 2 lnr 
The density of states thus takes the form 



(29) 



e |e u l')+e u W]/2 



-KlnLr 



1 



(30) 



of the partition function of a single ± vortex in a log- 
arithmically correlated random potential, known to be 
related to the one of a directed polymer on a Cayley tree 
0,0. A simple average of the partition function Zl, i.e. 
p e (0) indeed leads to Eq. (|27|), however, as is well known 
in the directed polymer problem only the logarithm of 
the partition function In Zl is self averaging. This im- 
mediately yields (2^j2^,|l^] our result for the typical DOS 
at finite size: 



Pt yp (0) ~ L 



-2+z 



(31) 



with z given by Eq. (poj). 

To identify the role of e we consider the first order 
terms e J d 2 r'(exp[±0(r) ± 0(r')). The typical value of 
each of these terms scales as the typical value of Zf/L 2 
(for opposite charges it is true in the massive phase we 
are interested in). This allows to identify a crossover 
function f(x), where 



Pt yP (Q) 



-2+z 



eL 



-2+2z 



L 



-2+z 



J{eL z ) (32) 



with f(x) = l+x at x — > 0. For i»l we can recover Eq. 
( p3| ) if the crossover function satisfies f(x) ~ (l/x) 1-2 ^ , 
i.e. the typical value pt yp (0) crosses over to the average 
p c (0) at e > L~ z , with both limits exhibiting a phase 
transition. 

This statistics can be described in a simple phenomeno- 
logical picture. A finite e provides a length scale (the 
vortex separation) and in effect cuts the system in inde- 
pendent pieces of sizes L e = e~ z . One has thus roughly: 



Pe(0) 



1 

L 2 



E ■ 



(33) 



(i) 

where the random variables Z L are independent with a 
lognormal distribution. For large (L/L e ), however, the 
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above sum acquires a normal distribution. A similar pic- 
ture was used to describe the related random diffusion 
problem, where the local first passage times are analo- 
gous to the local DOS in the present problem, and an 
external force produces a finite length scale. Analysis of 
the various regimes is described there and are expected 
to be quite similar here. 



C. relation to random diffusion models 

It is instructive to compare our results to the one ob- 
tained for random diffusion problems. As mentionned in 
the introduction, general random Dirac problems can be 
mapped onto random diffusion operators, which in gen- 
eral may involve two species, with absorption, creation 
and transformation. It is particularly simple in the case 
W = (random vector potential alone) where it maps 
onto a random Fokker Planck diffusion operator of the 
type: 



H FP P = V 2 P - V(F T + F L )P = -E'P 



(34) 



which describes the Langevin diffusion of a particle, 
dtP = HppP where P(r) is the probability that the 
particle is at point r at time t, in the presence of both a 
potential random force Fl = — Vt/ and a transverse one 
(a random drift), with divFx — 0. Equivalently, setting 
P = e- u / 3 if>: 

H' FP ip = (V 2 - F T V - (W) 2 + V 2 V)ij = -E'tp (35) 

with V — U/2 (K — 1 here). The operators Hpp and 
H' F p have the same spectrum. In two dimension, taking 
the square of (Q) with W = yields: 



-Hi 



V 2 -(Al+A 2 y ) 



a z (d y A x 



-2MA • V - iV -A 



(36) 
(37) 



with A 



1 identical to H' FP with A x 



d y V and 



-d x V (and V 



-V for the other component 



of the spinor) and Ft = 2iA. This is thus Arrhcnius 
diffusion in the random potential U with an additional 
imaginary random drift p5| . The diffusion dynamical 
exponent Zd associated with H' FP should thus be simply 
Zd = 2z. Note that all the operators obtained by vary- 
ing A have identical ground state wavefunction, ip ~ e~ v 
since the additional drift term does vanish in the ground 
state (in the diffusion context it means that the drift is 
along equipotentials of U). It is thus reasonable to expect 
that each of these models are described by a line of fixed 
points and that they all do exhibit a freezing transition 
for any value of A at the same value of a = ath = 2. 

In the absence of this additional drift (i.e. setting 
A = 0), the problem reduces to the one studied in J24] 
where indeed it was found that there is also a freezing 



transition in the dynamical exponent Zd in d = 1 and 
d = 2 with (assuming conventional dynamical scaling) : 

z d (X = 0) = 2 + 2{a/a th ) a < a th (38) 
z d (X = 0) = 4-v/ct/W a- > oth (39) 

Although it does indeed exhibit a freezing transition 
at the same threshold <j t h = d one sees from () that 
Zd(X = 0) < Zd(X = I) = 2z, i.e. the imaginary drift slow 
down the diffusion, presumably through interference ef- 
fects. It would be of interest to determine Zd(X) as well 
as to study freezing transitions in a generalized class of 
these diffusion models in two dimensions. 

It is possible to consider various one dimensional re- 
striction of the Dirac model, e.g. the so called super- 
symmetric quantum mechanics which also exhibits band 
center derealization p(i| , |3l"| | . With a log correlated U(x) 
this model was studied analytically in |24j , thus we know 
in that case the exact z = Zd/2 dynamical exponent of 
the random Dirac operator. 

IV. FULL QUANTUM HALL PROBLEM 
A. one layer problem: scaling of the DOS 

Finally, we consider the Dirac model where the scalar 
random potential in Eq. (^|) is retained, which describes 
the full quantum Hall system. We will first determine the 
DOS first at zero energy, and later around zero energy. 

In presence of a random scalar potential one has 

Y[n] - (iW - e)£ n « in Eqs. We can safely 

set e = in the definition p e (0) ~ dF var /de since W(r) 
provides a mass parameter. The variational method is 
similar to the previous case, except that e is replaced by 
iW(r) in Eq. (p"8|). Since the integral is dominated by 
large u and w (for ui > 0), it has the form: 



iWe 



U-\-UJ 



" (1 + iWe u +^) 2 



l (l + W 2 e 2u+2LJ ) 2 



(40) 
(41) 



using the ± symmetry of the W average |33|]. The lat- 
ter form is equivalent to the previous integral Eq. ( pc| ) 
with e replaced by the disorder average (W 2 ) — 8 (for 
the starting QH system 6 ~ a) and K is replaced by 2K . 
Hence a r ~ S 2 ^ z where now 



z =2-2K + 4aK z 
z' = 2K{V8a-l) 



a < 1/2K 2 
a > 1/2X 2 . 



(42) 



The DOS at zero energy can be written as 

pU+U 



p(E = 0) ~ (cos 6) ~ { 1 + w2e2u+2u )„,w ~ S a (43) 



where 



6 



2 

7' 



z 



(44) 



Since z' has a transition at a = 1/2K 2 the DOS has two 
transitions, at a = 1/2 and at a = 2 (for if = 1). The 
exponent a in Eq. (44) is the one expected from a scal- 
ing form p e ,,5(0) = 8 a g{e/8 z l z ) which connects with the 
6 = case solved in Section III (which requires g[x) — > 1 
at x — > and g(x) ~ x az / z at x — > oo). As we will see 
below there are however three phases, each with a differ- 
ent scaling function g{x). Note that z'/z increases from 
at small a to z'/z — 2 at a = 2 and remains equal to 
this value for stronger disorder. 

We have thus shown that the DOS is finite at the quan- 
tum Hall transition, which is a well established result . 
In addition however, we have determined the scaling of 
the zero energy density of states with the scalar poten- 
tial strength in the random Dirac problem. Note that 
in the previous case with e — > renormalization of K 
and a was of higher order in e and could be neglected. 
Here, however, a finite (W 2 (r)) renormalizes both K and 
cr, with a flowing to stronger values. This does not spoil 
our result for the exponent z' in the limit S — > however, 
as it would also yield only subleading corrections in 5. 

The physics of the problem becomes more apparent 
when one notes that the random scalar potential pro- 
duces (see e.g. Eq. (^) ) a random (imaginary) fugacity 
with a random sign. One can then again either extend the 
RG of Jl6| to this situation or consider the extreme dilute 
limit (single vortex) as in Eq. (l2gj). Both considerations 
lead immediately to a mapping onto the directed poly- 
mer on the Cayley tree with random Boltzmann weights 
of random sign (in the bare model the random sign acts 
only at the leafs of the tree but since both signs have 
equal probability one easily sees that this is equivalent to 
a random sign on each branch of the tree). This model 
was also solved in [34f| and is known to indeed exhibit a 
transition at half the value of the same sign problem, due 
to interference effects which in effect bind two replicas. 
The value of the z' exponent obtained by this method 
is identical to the one given above (Eq. J42])) . It is re- 
markable, and encouraging, that the variational method 
also captures this physics. 

We can now extend these considerations and obtain, in 
implicit form, the crossover function which describes the 
DOS p(E) in the small 5 limit. We will obtain explicitly 
p £ (0), from which p(E) can be extracted inverting 

The equation which determines a c and the DOS is now: 



dLU — — 7TTTK — T7T- ) OJ ,W 



1 + 2ee u+ul + (W 2 + e 2 )e 2u+2ul 1 



p e (0) = ( du 



, ;I 2ii+2^ 



1 + 2ee u+UJ + (W 2 



(45) 



(46) 



the first line is really o~ c +o~q but we can use that o~o is sub- 
dominant and only dominant exponents for the uj > are 
retained [p3|. We immediately see that there are several 



regimes, according to whether one can neglect all denom- 
inators (weak disorder), or whether the averages will be 
dominated by the rare events where either the terms pro- 
portional to W 2 or to e or both, are of order one. There 
are in fact three phases, in each of them scaling holds 
with different scaling functions f ,g,1Z: 



o c = 6 2 ' z ' f{e/5 z ' z ') 



(47) 
(48) 



from which one can extract the DOS scaling function 

p(E) =E 2/z - 1 n{E/5 z/z ') (49) 
determined implicitly by the relation: 



g(x) = - f dyy^ny)-^ 
7T J x A + y 1 



(50) 



weak disorder phase a < 1/2K 2 : 

Neglecting all denominators the equation for a c and 
p £ (0) become: 



pM = *l-*' 2 

l = eo-- z / 2 + 6v- z '/ 2 



(51) 
(52) 



which yield the scaling functions in implicit form: 

x = /(x)*/ 2 (l - /(x)-*'/2) (53) 
g{x) = f(xf- z l 2 (54) 

possibly the exact ones, up to prefactors. 
strong disorder phase I 1/2K 2 < a < 2K 2 : 
There one can still neglect denominators in averages 
involving the terms proportional to e but not in the one 
involving the terms proportional to W 2 . One finds: 



77TT 1— a/2 



Or = eer„ ' + a c 



(55) 
(56) 



Simple expansion shows that scaling still holds but the 
scaling functions are now implicitly given by: 



x = x = f(x) z ' 2 {l-f{x)- z / 2 ) 

*' = 4 (!--i=) 



g(x) = f{xY- z ' 2 



(57) 
(58) 

(59) 



This scaling function is accurate only up to loga- 
rithmic prefactors. A more accurate form is a c ~ 
5 2 / z '|ln5|- 2 T/ J 7(e^ z/z '|ln<5| 27/5 ') with 7 = 1/2 but 7 
is likely to be corrected upon a more careful RG treat- 
ment. 

strong disorder phase II 2/K 2 < a: 

At even stronger disorder, since z' = 2z, we expect the 
scaling region to be e 2 ~ S. To show that this is the case 
and to get some approximation for /(x) we notice that 
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the equation for a c can be approximated in the scaling 
region by: 

a c « Proba(ee u+w + (e 2 + 5)e 2{u+bj) > 1) (60) 

Solving the quadratic equation in e u+u yields that a c is 
indeed of the form (^), with some form for f(x), which 
here is approximate. We have not attempted to solve 
more precisely for f(x) or g(x) in this phase. 

Finally note that this crossover can also be studied at 
finite size, and is there complicated as it will probably be 
described as in p3 by a non trivial phase diagram. 



B. Transport in the quantum Hall system : the two 
layer problem 

We address now the more difficult problem of describ- 
ing transport and localization in the QH system. Trans- 
port is derived by a disorder average of advanced and 
retarded propagators, hence two partition sums corre- 
sponding to the action of Eq. (||) with ±e. The problem 
is then of two layers with fields 6-\ and 0i and common 
disorder A(r) and W(r). The role of ±e is to deter- 
mine the proper ground state near which the variational 
method applies, i.e. for — e we shift 9[ — > 6^ + it so that 
the nonlinear terms become 

exp{-^[(-W + e)cos0 T)O + (iW + e) cos0j. )O ]} (61) 

a 

where a — 1, m are replica indices for each layer and 
we have redefined here W/ira, e/ira — > W, e, respectively. 
Expansion in powers of iW ± e and identifying the dom- 
inant terms of the form exp(in • 9) yields 

Him = ^(-iW)£a n ^iW)£° n *°exp(in-0) (62) 



with e now set to zero, n is now a vector of 2m entries 
(n-fi, ...n-\ m , nil, nj, m ) and similarly for 0. 

The Gaussian part of the Hamiltonian can be written 

as 

W' = £^(|0 + , a (q)| 2 + |0-,a(q)| 2 ) 



q , a 



E ^^ + ,a(q)^,6(-q) 

q,a,6 



(63) 



where 9± a = ($^ a ± 6j, )/v2- Since only the 0+ mode is 
affected by the common disorder A(r) we expect 9± to 
have distinct self masses. The variational Hamiltonian is 
then 

H = H' + ^ [<rt\ e ±a(<l)\ 2 5 ab + a£6 ±a (q)6 ±b (-q)] 

q,±,a,6 

(64) 



The propagators of the ± modes are used to define 

5^<0±a(q)0±i(-q)>o = -^u ± 8 ab - 2A± 
q 

u ± = -{K/A) \n(A c /4:TTKat) 

A+ = <jK 2 ln(A c /47r!fcr+) + Ka+/2a+ - aK 2 

A- = Kao/2*-. (65) 

We write the interaction Eq. (|62l) in the form 

(H int ) = ((-iwj£»<(iwj£« n °° 
exp[u+ ^(n Ta + n ia ) 2 + u + J^(n Ta + n [a ) 

a a 

+U- E( n Ta - n iaf + U- E("Ta ~ n ia)])u (66) 
a a 

where the u> average reproduces the required form with 

^±(Ea n Ta±"-la) 2 , 



exp[ 



2A+ 2A- 1 2tt^JA+A- 



(67) 



The sum in Eq. (pq ) can be written as (Hmt)o = (Z m ) L 
with 

+W 2 e 4u++2ul+ + W 2 e iu - +2u; - 



+W 2 e iu -~ 2uJ - +W 2 e 4u+ ~ 2uJ+ 



(68) 



The masses are to be found by minimizing the variational 
free energy F var = F + (H — H )o, as in section III A. 
We expect to find a massless solution, e.g. crT = 0. This 
is indeed a possible solution with u_ — > — oo and 



Z = 1 + W 2 e 4u++2uJ+ + w 2 e 4u+ ~ 2uJ+ 



(69) 



VF 2 can be replaced by its average and then this has the 
same structure as the single layer problem of section IIIA 
with the replacement K — > 2K, i.e. the phase diagram is 
the same as Fig. 1 with the 1/K axis replaced by 1/2K. 
The starting line K = 1 is now tangent to the phase 
boundary at a = 0. However, for a > 9+ is massive. 

To find the QH localization exponent we introduce the 
mass term which corresponds to 



A d (sin(9 T -sinflj) = 2A d cos(<9+/V2) sin(0_/V2) (70) 

Note the opposite signs due to the shift of 9i as required 
by the sign of e. Since 9 + is massive, can be replaced 
by A d = 2A d (cos(6>+/%/2)) and defining 9 = #_/V2+tt/2 
leads to an effective Hamiltonian with 

H ef f = / d M^|V ■ 9(r)\ 2 A d cos0(r)] (71) 
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This system has a correlation length £, related to a mass 
£-2 ~ A^ which from first order RG is v = 2/(4 - if). 
For the starting QH problem with K — 1 the localization 
exponent is then v — 2/3 (the numerically known value 
is w 2.3). We expect, however, that iT is RG driven to 
a different value. This is beyond the variational scheme 
which gives reliable exponents only for small (W, e) cou- 
plings (but arbitrary a, K), as in the e — > case. E.g., to 
2nd order in replicated sine Gordon RG, the most rele- 
vant operator W cos(9 a + 8b) yields 

— = (2-2K + AaK 2 )W 
di 

^£l = -2{K-2aK 2 )W 2 

^ = (K - 2<tK 2 )W 2 (72) 

which shows that indeed K increases at weak disorder. 

We have also looked for other solutions of the varia- 
tional scheme of the form a~ ~ (o'^) a and found that 
only a = 1 exists. This corresponds to decoupled layers 
with the phase diagram of Fig. 1 (K replaced by 2K). 
Hence at K = 1 this is a massive phase and does not 
correspond to the QH problem. 

Our method shows that a subset of the degrees of free- 
dom can form a massless phase and we hope that it will 
stimulate further progress. 

V. CONCLUSION 

In conclusion, we have derived an exact formulation 
of two dimensional random Dirac fermions in terms of a 
random sine Gordon model, or equivalently a disordered 
Coulomb gas. The DOS of the Dirac fermion system 
identifies with the expectation value of a cosine operator 
in the sine Gordon model ( equivalently the charge fu- 
gacity in the Coulomb gas ) and the dynamical exponent 
z as its scaling dimension. We found that at zero energy 
and with random vector potential only the Dirac system 
maps onto the random gauge XY model with infinitesi- 
mal vortex fugacity (dilute limit of the CG). Using meth- 
ods and results from previous studies of the XY model, 
we have computed the exact dynamical exponent z for 
the random vector potential model at any disorder, and 
thus obtained the critical behavior of the DOS around 
zero energy. We found that it exhibits a transition at the 
same threshold value than the previously known transi- 
tion in an exact ground state wave function. It shows 
that all eigenstates near the band center are affected by 
this transition. The physics of this transition is closely 
related to the freezing transition in the XY model in the 
limit where the vortex core energy is taken very large. As 
we show here the density of states in a finite size sample 
becomes broadly distributed with typical values scaling 



differently than average ones with the system size. This 
corresponds to the eigenstates being peaked around some 
few centers in the sample. 

It is likely that similar freezing phenomena are of im- 
portance in a broader class of two dimensional disordered 
models. They were recently found to occur in random 
diffusion models, for instance in the problem of random 
Arrhenius diffusion of a particle in a logarithmically cor- 
related potential. We have analyzed the similarities and 
differences of the transitions which occur in the dynam- 
ical exponents z of both models. As was shown in |24| ] 
strong disorder renormalization group captures the dy- 
namical behaviour in the glass phase, which suggests that 
it could be used to study the present problem as well. 

It is of utmost importance to understand what are the 
consequences of this transition when a random potential 
is added, corresponding to the quantum Hall system. As 
a first step in that direction we have determined the scal- 
ing of the DOS (which is then finite at zero energy) when 
the additional random potential is small. There we found 
two transitions, one at the same value than the random 
vector potential model, the other at a much smaller value 
of disorder. It would be nice to know whether these tran- 
sition lines extend away from the random vector potential 
fixed line. Numerical checks of our predictions as well as 
a numerical calculation of a glass order parameter (e.g. 
Hx IVK^OI 4 ) m the full model would help understand 
these issues. 

We believe that such freezing phenomena, originally 
studied in the context of disordered Coulomb gas and 
XY models, will also affect a broader class of disordered 
fermion models in two dimensions. This can be studied 
systematically by extending the bosonization approach 
introduced in the present paper, and for instance, search- 
ing for all perturbations around the random gauge fixed 
plane and computing their (non trivial) scaling dimension 
(here it was done only for the vortex fugacity operator). 
In particular it is of interest to know whether the non 
linear sigma models studied in || also exhibit freezing 
phenomena. They are indeed generalizations of the Li- 
ouvillc model which captures the single vortex problem 
and does exhibit a freezing transition. Results are al- 
ready known in related cases. For instance it was shown 
in flU that if one adds pinning disorder to the random 
gauge XY model the vortex density (DOS) acquires a 
~ exp(— | lne| 2 / 3 ) dependence in the vortex core energy 
ln(l/e) with a non trivial exponent. 

Finally, we have formulated the problem of the trans- 
port in the quantum Hall system as two coupled random 
sine Gordon models. We have applied the variational 
method, which should capture some of the non pertur- 
bative effects. We have analyzed the system using in the 
plane of two (dimensionless) parameters and were able 
to find a massless phase. Though qualitatively encour- 
aging, to obtain the quantitative characteristics of this 
phase requires to incorporate in a more precise way ad- 
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ditional renormalizations of these two parameters. It is 
tantalizing that possible values of these renormalized pa- 
rameters (consistent with numerics) seem to lie within 
the region near the glass phase boundary. 

We thank Bernard Derrida and Leon Balents for use- 
ful discussions. B.H. acknowledges support by THE IS- 
RAEL SCIENCE FOUNDATION founded by the Israel 
Academy of Science and Humanities. 
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